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ABSTRACT 



A model for calculating the index of refraction of 
atomic species present in the cavity of a Xenon Fluoride 
laser is applied to Neon and Neon* (3s). The model con- 
siders the variation from unity to be a function of the 
absorption cross section. Below the ionization threshold, 
the cross section is a set of discrete transitions between 
the various energy levels of the specie of concern. Above 
the 5p level, the discrete transitions become indistinct 
and the continmim cross section function is extended into 
this region. The continuum cross section function is a 
polynomial function fitted to independent research results. 
The results for Neon I agree with previous theoretical and 
experimental results. The model is extended to Neon* (3s). 
Comparison with the results of independent research of the 
polarizability and cavity phase shift derived from the 
calculation indicate accuracy is within ten percent. 
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I. 



INTRODUCTION 



The earliest successful high powered lasers operated 
in the vibration/rotation transitions of CO 2 and produced 
radiation in the infrared, principally in the 10.6 micron 
line. Extensive research and development of technology in 
the infrared laser field has produced impressive results 
that clearly demonstrate the feasibility of directed energy 
weapons utilizing these transitions for laser power produc- 
tion. The relatively long wavelength results in large 
optical components to form and direct the beam. Even though 
the highest power lasers are of the CO 2 type, their ex- 
tremely large size limits their usefulness to laboratory 
experiments. 

The shorter wavelength of chemical DF lasers provided 
improved efficiency and reduction in size of the power 
production and cavity components. This was the first laser 
system to demonstrate target engagement and destruction.^ 
Still governed by the laws of diffraction, the optics were 
large and cumbersome. 

Early work on the excimer lasers had demonstrated their 
ability to lase in the vacuum ultraviolet; however, problems 
with mirrors and competing absorption processes lead to low 
efficiencies for these lasers. The channeling of atomic 



1. "International Defense Digest," International Defense 
Review, 12:316, Vol. 3, 1979... 
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ionization and excitation into these excited states was 

found to be very efficient. Further research by Velasco and 

Setzer produced interesting fluorescence in rare-gas halide 

combinations that led Ewing and Brau to the demonstration of 

the rare-gas halide exciplex laser, which is now included in 

2 

the broad "excimer’' description. 

The rare-gas monohalides are simple diatomic molecules 
which are bound in the excited state but are unstable or 
slightly bound in the ground state. The excited state of 
the excimer species can radiate in a broad band that is red 
shifted from the wavelength of the parent excitation. These 
lasers have demonstrated electrical efficiency as high as 
ten percent in various wavelengths, as depicted in Table I. 



TABLE 1. Rare-gas Halide Laser Emission Wavelengths 
Peak wavelengths of most intense band (nm) 





Ne 


Ar 


Kr 


Xe 


F 


107^ 


193 


248 


351,353 


Cl 


b 


170 


222 


308 


Br 


b 


166^’^ 


206^ 


282 


I 


b 


b 


185^’^ 


252^ 



3 . 

( Has not demonstrated laser oscillation.) 

( ^Predissociates , hence emission is weak or unobservable.) 



2. Ewing, J. J. , Rare-Gas Halide Lasers , Physics Today, 
May 1978, pp. 32-39. 
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Attempts to realize the inherent advantages of high 
electrical efficiency at short wavelength have resulted in 
extensive research and development activity in this area. 
Extensive literature on the kinematics of the excitation 



There is a demonstrated understanding of the chemical kin- 
etic processes involved in energy transfer and formation 
of the excited lasing molecules. Engineering development 
has demonstrated improved flow and acoustics effects. The 

latter were essential to the development of good beam quality 

0 

in a demonstration device. In this case, good beam quality 
is considered to mean uniform phase across the cavity aper- 
ture. Nonuniform phase causes the peak axis intensity to 
fall off according to 



shifts across the aperture. The spatial phase relation is 



3. Huestis, D. L., Small Scale Laser and Kinetics 
Program , paper presented at Rocketdyne International, Canoga 
Park, CA, 15 February 1979. 

4. Hogge, H. D. and Crow, S. C., Flow and Acoustics in 
Pulsed Excimer Lasers , paper presented at AIAA Conference on 
Fluid Dynamics of High Power Lasers, I 1-4, Cambridge, MA, 

2 November 1979. 

5. Hogge, H. Dwight and Crow, Stephen, Flow Design 
Concepts for Pulsed Visible Wavelength Lasers, Darpa Report, 
May, 1978. 

6. Siegler, R. , et al . , XeF Demonstrator Program , 
Rocketdyne International, 15 February 1979. 



3 4 5 

process and associated flow effects is available. ’ ’ 



= exp 
o 



where <t> is obtained by squaring and averaging the phase 



rms 
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ELECTRON BEAM 

PUMP 




The initial efforts to limit phase distortion considered the 
index of refraction in terms of the Gladstone-Dale formula- 
tion 

n = 1 + Gp 

where G is the Gladstone-Dale constant, hence the concern 

7 

for flow and acoustic effects. Hogge and Crow state that 

flow uniformity and reduced acoustic effects have achieved 
-5 

Ap/p = 5x10 , the required flow uniformity. 

It has been suggested, however, that the requirement for 
flow uniformity should be stated in terms of An/n due to 

g 

nonuniform pumping in the cavity. In a typical flow mixture. 



7. Hogge, ibid. 

8. Fuhs, A., Etchechury, J. , and Cole, L. , Progress 
Report # 2 , XeF Beam Quality Considerations , submitted to 
Dr. J. A. Mangano, DARPA, 25 February 1979. 
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the constituents of the medium (NerXetNF^) are present in a 

9 

fixed ratio (750:3:1). Even in a gas with fixed chemical 
composition, pressure and translational temperature, and 
hence, fixed Ap/p, An/n can vary according to the ionic and 
electronic excitation of the medi\im caused by pumping. 

Consider the Ladenburg formulation of the Krame^s-Kronig 
dispersxon relation 



n-1 = 



e^N 

Sirm^c-'^ 



N 



S/SI 






kj 



■kj 



N, 



1 - 






I-l 



There are t distinct species indicating specific electronic 
or ionic excitation present in the cavity. Index j sums 
over lower levels and k over upper levels of possible con- 
figurations. Identify the mole fraction as 



X 



£ ■ 




1-2 



the fraction of the total molecular population represented 
by species £. Finally, define a quantity 



2 x" 

T. = -g- V V 

^ 2im c^ ^ X^ - X 

e i k ^ 



2 ^kj 
kj 



1 

""j ^k 



1-3 



In terms of the quantities defined above 

n - 1 = 

I 



1-4 



9. Siegler, ibid. 

10. Mitchell, A. C. G. , and M. W. Zemansky, Resonance 
Radiation and Excited Atoms, p. 238, Cambridge, 1934. 
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or the index of refraction in the cavity is the sum of the 
contributions of all the species present in the cavity. 
Since = NX^, 

n - 1 = + NgT 2 + ... 1-5 



The partial density is given by = 
the molecular weight of tth species^ 
density, equation (1-5) becomes 

h - 1 = 



N^ra^, where m^ is 
In terms of partial 



+ P 2 '^ 2^^2 '*’■ ■ ‘ 



Introduce the mass fraction for the gas mixture 

= p^/p, where p is the density of the gas mixture. 

The result is 

h - 1 = p[Y^T^/m^ + +...] 1-7 

Equation (1-7) shows that n - 1 is proportional to p; 
however, equation (1-7) further demonstrates that n - 1 
depends on mass fraction and index of refraction of each 
species. Even with uniform flow, the spatial phase relation 
across the cavity will not be uniform. Accurate modeling 
of the phase effects requires knowledge of the populations 
and the contribution to the index of refraction of the 
species in the cavity. This paper will explore a model 
for the index of refraction of the atomic species in the 
cavity and apply the model to neon. 

In Section II, the formulation of the index of refrac- 
tion to be applied in the current work will be derived. 

In Section HI, the concept of transition arrays and the 
calculations of the transition array oscillator strengths 
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will be reviewed. The available data for the cross section 
of Neon I (2p®) and Neon* (3s) will be reviewed in Section 
IV. The formulation derived in Section II will be applied 
to Neon I and the results compared to previous theoretical 
and experimental work in Section V. Section VI will discuss 
the calculations for Neon*(3s) and these results will be 
compared to independent experimental results in Section VII. 
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II. INDEX OF REFRACTION FROM DISPERSION THEORY 

The ratio of the speed of light in a vacuum to the 
speed of light in a medium is known as the index of re- 
fraction, 

c 

n = V II-l 

where v is the phase velocity. Consider a simple wave 



where the phase relation is determined by 

03t-kx = constant 



Then 



V 



dx 

dt 




Substitution of (4) into (1) yields 



kc 

(0 



II-2 



II-3 



II-4 



II-5 



However, the index of refraction and propagation vector 
are frequency dependent, hence 

nC») = . 1I-5A 

^ (0 



In terms of the relative permittivity and relative 
permeability of the medium 



/£_ . y. 



/ 






For gas lasers, as for most non-ferromagnetic 
does not vary significantly from unity. 



II-6 

elements , 
This fact 
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iTiiiWi’' 



yields Maxwell's Relation 



n = /K = /K (co) II-7 

e e 

where is the static dielectric constant of the 
medium at frequency oj. 

The physical basis for the frequency dependence of n 
and the mechanism whereby the phase velocity in the 
medium varies slightly from c lies in the interaction of 
an incident electro-magnetic wave with the array of atoms 
constituting the dielectric material. In a gas, the 
molecules are sufficiently far apart that each molecule 
may be considered independent. The external field sep- 
arates positive and negative charges which creates an 
additional field component known as the electric polar- 
ization (dipole moment per unit volume) 

P = Nex 1 1-8 

where N is the number of molecules per unit volume and 
X is the average relative displacement of charge due 
to the external field. The following relations from 
Maxwell's equations hold 

D = eE = £ K,:.E = £ E + P 
o ^ o 

Simply stated, the electric displacement is proportional 
to the external field. We can write 

Ke = 1 ^ n-10 

5 = /k; = 

o 

11. Hecht, E. and Zajac, A., Optics , Addison-Wesley 
Publishing Company, 1974. 
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where n is the complex index of refraction. 



n^ = 1 + ^ = 1 + N ex 



11-12 



£ E 
o 



£ E 
o 



For more than one dipole (as in array of N atoms) 



5 = S IV 11-13 

N ° 

or the total change in the index of refraction is equal 
to the sum of the contributions of the dipoles. 

An analytical expression for the average displacement 
of the charges is easily developed in classical terms of 
the Lorentz oscillator model of the atom. Although dev- 
elopment of this model is properly done only with quantum 
mechanics, the classical treatment gives similar results. 
In the Lorentz model, the outer electrons are considered 

bound to a nuclear core by a restoring force proportional 

2 

to the displacement from equilibrium x)* At 

equilibr ixim, the equation of motion for an electron is 



“e ^ “e’f af ° 

dt 

where is the natural resonance of the system. If 

the electrons are acted upon by a harmonic electromagnetic 

io)t 

wave of the form E = E^e , the system resembles a 
classical forced oscillator and Newton's second law yields 



d X , dx , 2 _ i(ot -iir 

e ITT + “e at ^ ^^'35 

dt 
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Solutions to equation (11-15) are easily found by assuming 
a solution of the form of the forcing function. They yield 



X 



e/nig 
2 ^2 

0) + oj + iyco 

o ' 




11-16 



Substituting into (11-13), and canceling the field terms, 
yields 

-2 V - 

"- 1^1 —2 2 — 

£ 03^ + oj + lym 



and we have the form of the complex index of refraction. 
Separating equation (11-17) into its real and imaginary 
components 

n = n + i/C 11-18 



gives 



n(m) = 1 + I X 
£ 



(N^e^/me^)(aj^2_ 

7 2 2T“7 2 2“ 

(m^ - 03 ) + 03 Y» 



,,, 1 ^ (N^e^/mc,^)(co Y^3 

K(o3) ~ 2 2- , 2 2. 2 



(03^ + 03 ) + 03 



11-19 



Kramers and Kronig first derived comparable equations 

from quantum mechanical considerations. Their equations 

are found in various texts, and the form published by 

12 

Landau is shown for comparison. 



12. Landau, L. D. , and Lifshitz, E. M. , Electro- 
dynamics of Continuous Media , translated from Russian by 
J. B. Sykes and J. J. Bell, p. 259, Addison-Wesley Pub- 
lishing Company, 1960. 
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n (( 0 ) - 1 + 






K(a)') dw'^- 

,'2 2 

0) ' - CO 



11-20 



K(oi) = - 



ip/- -1 doj' 

" y o ,''2 - 0.2 



11-21 



where P stands for the Cauchy principle value. Using 
the fact that K(co) is an odd function 



n(co) = 1 + 



1 / 



co‘ /C(o)) dco" 

,2 2 

CO - CO 



0 



(11-18) can be rewritten in a more convenient form 



11-22 



where 



or 



and 



6(E) 



n(co) = n(co) + ifC(co) 



n(co) = 1 + 5(co) 



n(E) = 1 + 5(E) 



^ 2 r E' K(E ' )dE 
^0 



II-23 



11-24 



11-25 



11-26 



where the equation is now considered in terms of the 
energy of the incident photons, E = "feco, and 6(E) is 
the variation of the index of refraction from unity. 
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13 

Levinger demonstrates the absorption coefficient K 
is related to the total cross section of the individual 
atoms by 



K(E) 



Ncfe 

2E 



cr(E) 



11-27 



or 



2E (E) = Nc*ha(E). 



II-27a 



Substituting into (11-26) 



5(E) 



Ncfi r” g(E)dE’ 
^ E'2 _ 



11-28 



and we now have an expression for the variation from 

unity of the real part of the index of refraction in 

terms of the total cross section of the atom and the 

14 

energy of the incident photon. Cooper states that the 
majority contribution to the index of refraction comes 
from photon energies significantly below 100 Rydberg 
(one Rydberg is 13.595 ev), and that in the region from 
zero to 100 Rydberg, the scattering cross section is 
negligible compared to the absorption cross section. 

Up to the ionization threshold, the absorption cross 
section is composed of a series of discrete transitions 
defined by the transition rules of quantum mechanics. 
Beyond the ionization threshold, the absorption cross 
section is a continuum to the limits of integration. For 



13. Levinger, J. S. , letter dated 4 April 1979 to 
James Etchechury. 

14. Cooper, J. W. , "Photoionization from Outer Atomic 
Subshell. A Model Study," Physical Review , 128:681, 15 
October 1962. 
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a discrete transition, the total absorption cross section is 
related to the oscillator strength for the line by 



y'a(E)dE = 

line 




me 



f 

nm 



11-29 



The oscillator density for a line is related to the 
cross section by 



df 



nm 



dE 



me 

O 2 2j- 

2it e h 



0(E) 



11-30 



wh:ch now defines the constant R 
rewrite (11-28) 



— 2 ■ 2 _ > 
2-n^e'^h 



where 



5(E) 



C 



^E ^-E^ '' E’2- E^ 



11-31 



Nch 

ttR 



E is the discrete transition energy 
nm 

E' is the variable of integration 
and- energy units are Rydbergs. 

15 

This formulation was used- by Migneron to calculate 

IS 

5(E) for Helium and by Liggett. to .calculate 5(E) for 



15. Migneron, R. and J. S. Levinger, "Index of 
Refraction and Sum Rules for Helium," Physical Review , 
13:A646-48, 2 August 1964. 

16.. Liggett, G. and. J. S. Levinger, "Calculation of 
the Index of Refraction of Neon and Argon," Journal of 
the Optical Society of America , 58:1:109-113, January 
1968. 
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Neon and Argon. Their theoretical results agree well with 
existing experiments, indicating good accuracy for this 
technique . 
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III. REVIEW OF THE OSCILLATOR STRENGTH 



In order to evaluate the formula developed for the index 
of refraction, the oscillator strengths in absorption for the 
respective transition arrays are necessary. Condon and 
Shortley define a transition array as the totality of lines 
resulting from transitions between configurations (e.g., 
ls^2s^2p^to ls^2s^2p^3s ) . Bethe and Salpeter^^ give the 
oscillator strength in a central field model for a single 
electron transition between configurations n and m, defined 
by quantum numbers nt and n' V . 



f 



1 

nm 5" 



(E -E )R ^ 
ZL+l n m nra 



For £’=£-1, £ =£, and for £'=£+1, £ =£+l 

’ max ’ ' max 

E^ are the energy eigenvalues of the states 

is the radial matrix element 
nm ~ 



, E and 
’ n 

n and m. 



III-l 



''nm =7 P^(r)dr III-2 

0 

The single electron oscillator must be multiplied by the 
number of electrons in the outermost shell to determine 
the absolute oscillator strength. 

Gruzdev has calculated the radial transition integrals 
for neon using a standard program which calculates the 



17. Cited in Cooper, J. W. , "Photoionization from Outer 
Atomic Subshells. A Model Study," Physical Review , 128:681, 
October 1962. 

18. Gruzdev, P. F., and A. V. Loginov, "Neon, Radiative 
Lifetimes for the Levels of 2p^ms, 2p^np, 2p^nd Cm = 3-6, 

n = 3-5) and 2p54f Configurations," Optics and Spectroscopy , 

35:1-3, July, 1973. 



Hartree-Fock radial wave equation for the configurations 
considered. These are reproduced in Table 2. 



TABLE 2: Transition Integrals / P^(r )rPj^(r )dr for Neon Atoms 



Transition 


HF 


Transition 


HF 


3s-2p 


0.43 


3d-2p 


0.11 


3s-3p 


4.76 


3d-3p 


6.3 


3s-4p 


0.35 


3d-4p 


9.22 


3s-5p 


0.13 


3d-5p 


0. 11 


4s-2p 


0.17 


4d-2p 


0. 081 


4s-3p 


4.25 


4d-3p 


1.6 


4s-4p 


10.91 


4d-4p 


11.4 


4s-5p 


1.07 


4d-5p 


19.2 


5s-2p 


0.10 


5d-2p 


0.06 


5s-3p 


0.98 


5d-3p 


0.82 


5s -4 p 


10.3 


5d-4p 


3.22 ' 


5s-5p 


19 .34 


5d-5p 


17.68 


6s-2p 


0.069 


4f-3d 


10.2 


6s-3p 


0.51 


4f-4d 


15.9 


6s -4p 


12.06 


4f-5d 


1.74 


6s-5p 


18.2 







For equation III-l, an average level found by summing 
over the levels of the configuration given in Moore will be 
used in place of the energy eigenvalues used by Cooper, i.e., 



where i indicates the various energy levels of the 
configuration and g^ reflects the degeneracy or quantxim 
statistical weight of the level. 
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While the energy eigenvalues could be found by evaluating 
the radial Schroedinger equation, Ligget's discussion of the 
calculation of the index of refraction cites difficulty in 
calculating energy eigenvalues. Use of published informa- 
tion for was considered to be more accurate and, in 
addition, expedited the calculations. 
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IV. THE ABSORPTION CROSS SECTION 



Extensive research data are available for Neon 
2 2 6 

(Is 2s 2p ). Since the demonstration of the Helium-Neon 

laser, a vast accumulation of theoretical and experimental 

work has occurred. It is safe to say that the important 

characteristics of the atom in its ground state are well 

known to researchers in this area. Not so well known are 

the characteristics of the metastable (3s) configuration; 

as a consequence, the accuracy of the available data for 

the photoionization cross section is not known. 

19 

Ederer and Tomboulian measured the total cross 
section of ground state neon from the ionization limit, 

1.58 Rydberg. Their results are tabulated in the reference 
and plotted in Figs. 2 and 3. In his experiment, Ederer 
measured the transmission of ultraviolet light with the 
gas present and with vacuum. The cross section is given 

*’5' , I 

where N is the density of atoms and L is the path 
length. To get monochromatic radiation, he used a normal 
incidence spectrograph. To measure the intensity of the 
current, both a photographic plate and a Geiger -Muller 
counter were used. 



19. Ederer, D. L. , and Tomboulian, D. H. , "Photoioniza- 
tion Cross Section of Neon in 80-600A Region,'^ Physical 
Review, 133 :6A: 1528, 16 March 1964. 
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CROSS SECTION OF Ne (2P6) 
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CROSS SECTION OF Ne‘ (2P6) 
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Figure 3. Photoionization Cross Section of Neon (2p 
in the Region 3.3 to 11 Rydberg. 



20 

For energies above 12 Rydberg, Piech and Levinger 
have fitted the cross section into a series of power laws. 
The power laws are fit to experimental data developed by 
other researchers. The laws and their respective limits 
are listed below: in Table 3: 



TABLE 3; Empirical Fit for 
Using Power Laws. 


Photoionization Cross Sections 


Energy (Ry) 


Cross Section (Mb) 


13.0-63.7 


1.12xl0^E"^-^^ 


63.7-227.8 


3.24xlO^E"^*^^ 


227.8-911.3 


1. lOxlO^E”^*^® 



Theoretical and experimental knowledge of the photo- 
ionization of metastable rar e-gas atoms is rather scarce. 

21 

McCann and Flannery have employed the single particle 
central field approximation to examine the photoionization 
of Ne*(3s). Their results are plotted in Fig. 4. The 
cross section of the outer s electrons of the metastable 
gas as a function of the ejection energy has been found to 
yield considerable improvement over standard calculations 
based on an interaction obtained from iteration of the 
ground state configuration. These calculations demonstrate 



20. Piech, K. R. , and Levinger, J. S., "Sum Rules for 
Neon Photoeffect," Physical Review , 135:A332-6, March 1964. 

21. McCann, K. J. , and Flannery, M. R. , "Photoioniza- 
tion of Metastable Rare-Gas Atoms," Applied Physics Letters , 
31:9:601, 1 November 1977. 
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that the cross section is much smaller than that of Neon I 



and present an optical window in the 2000-3000 A region. 
Although not in the immediate region of consideration 
(3510-3530 A), the window is close enough to have consider- 
able effect. Further inner shell ionization occurs only 
when relatively short wavelengths (@ 600 A) are reached. 
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CROSS SECTION Ne» (3s) 




(2-20 I *aw) N0I1D2S SSOaD 



32 



Figure 4. Photoionization Cross Section of Metastable Neon (3s) 



V. CALCULATIONS OF THE INDEX OF 
REFRACTION OF NEON (2p6) 

Equation (11-31) was used to evaluate the index of 

refraction in the region of interest (3510-3530 A, or 

.258 Ry), at standard temperature and pressure. If the 

—18 2 

cross section is measured in Megabarns (1 Mb=10~ cm ) 
and the energies are considered in Rydbergs (1 Rydberg= 
13.595 eV), then the respective value of the constants 
in equation (11-31) are: 

C = 1.00009x10"^ Ry^ 

R = .12895 Ry"^ Mb"^ 

The oscillator strengths for the 2p->(3s, 4s, 5s, 3d) 
transitions were calculated using Gruzdev's transition 
integrals. For the region of interest, the contribution 
to the index of refraction by inner shell ionization and 
any double excitation is negligible. The calculations 
agree well with the calculations of Cooper, except for 
the 2p->3d transition. However, this discrepancy between 
Cooper and Gruzdev does not adversely affect the results 
of the calculation due to the relatively small influence 
the 2p^-3d transition has on the oscillator sum term of 
equation (11-31). The results of the calculations are 
given in Table 4. 

The discrete cross sections above the 5s or 3d level 
become merged and insignificant relative to the numerical 
values of the lower transitions. This is clearly demon- 
strated in the calculations of Table 4 and in Figure 5. 
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Region I Region II 



Figure 
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AH known lev eh included 
in the system 
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5. Energy Levels of Neon I and the Region 
of Discrete Oscillators Arrays (Region 
I) and Region of Extension of the 
Continuum Integral (Region II). (Scale 
in thousands of cm~l.) (From Condon 
and Short ley, p. 302.) 
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Table 4. 


Tr ansition 
for Ground 


Array 

State 


Oscillator 
Neon (2p®) 


Strengths 


Transition 


Etchechury 




Cooper 


Process 


2p->3s 


0.1632 




0.1632 


Absorption 


4s 


0.0279 




0.0264 


Absorption 


5s 


0. 0102 




0.0090 


Absorption 


3d 


0.0239 




0.0372 


Absorption 



To account for transitions to 5s and above, the continuum 
cross section data were extrapolated into this region; the 
limit of integration was adjusted to 1.51 Ry, which is 
determined by energy of the 5s level. 

Adjusting the limit of integration in this manner is 
considered valid, since the discrete oscillator strength 
distribution can be considered as an extension of the con- 
tinuum oscillator strength distribution. Cooper discusses 

23 

this concept in greater detail. 

In the continuum, a polynomial was fitted to Ederer's 
experimental cross section data in the region 1.58 to 3.3 
Ry and a separate polynomial function from 3.3 to 13.0 Ry. 

This resulted in a set of polynomial functions that, when 
separated into individual components, resembled the form of 
the power laws developed by Piech for the region above 13.0 Ry. 



23. Cooper, J. W. , ibid. 
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The integral could be represented as a set of separate terms 



/ 



g(E)dE* 

E'^-E^ 



th 



3.3 ,13.0 63.7 227.8 911.8 

^ f * f \f * I *f 

1.51 3.3 13.0 63.7 227.8 



where the integral has been terminated at the upper limit 
of 911.8 Ry. 

From Figs. 1 and 2 and Table 4, it is apparent that the 
major contribution to the index of refraction is made by the 
region near the ionization threshold. In the case of Neon, 
the contribution to the index from the region 1.51 Ry to 
3.3 Ry, characterized by the peak in the cross section curve 
is seventy-three percent of the continuum term. The region 
3.3 Ry to 13.0 Ry contributes twenty-seven percent and the 
contribution from the region above this is so insignificant 
as to not be considered in these calculations. This does 
not result in any significant error and demonstrates that 
the majority of the photoionization to the index of refrac- 
tion originates in the region near the ionization threshold. 
The small contribution at high energy is illustrated by the 
exceedingly small cross section that can be presumed by 
extrapolation of Fig. 3 to high energy. Comparison of 
results in Table 5 indicates good agreement for the 
Etchechury model, the more exhaustive calculations of 
Liggett, and the definitive experimental work of Cuthbertson 
The results demonstrate that the simplification of the 
problem results in small error and the accuracy of the 
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Etchechury model is considered sufficient for application 
to the electronic species Ne*. 





Table 5: Calculation of 


6(E) = 


= n-1 






Author 


Source 


(n-1 )xlO"^ 


E 




X 




Etchechury 


Theoretical 


0.6705 


0.258 


Ry 


3520 


O 

A 


Liggett 


Theoretical 


0.6778 


0.25 


Ry 


3630 


A 




Theoretical 


0.6825 


0.30 


Ry 


3020 


A 


Cuthbertson 


Experimental 


0.6788 


0.25 


Ry 


3664 


A 




Experimental 


0.6812 


0.27 


Ry 


3342 


A 
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VI. CALCULATION OF THE INDEX OF REFRACTION OF 
NEON* (3s) 

Equation (11-31) was used to evaluate the index of 
refraction of metastable Ne*(3s) at .258 Ry (3520 A) and 
typical XeF cavity conditions. The conditions were standard 
temperature and pressure and an electron temperature of 
13500°K due to pumping and kinetic action in the cavity. 

Using Boltzmann statistics with the assumed electron tem- 
perature of 13500°K, the ratio N(Ne* ) /N(NeI ) is 1x10"^. 

The respective values of the constants in equation (11-31) 
are: 

C = 1.00009 X 10"^° Ry^ 

R = 0. 12895 Ry~^ Mb”^ 

Drawing on the previous section, the Etchechury model 

was used to evaluate the index of refraction of Ne*(3s) in 

the region of interest at standard temperature and pressure. 

The oscillator strengths were calculated from Gruzdev's 

integrals for the 3s->-(2p, 3p, 4p, 5p) transitions, as listed 

24 

in Table 6. McCann and Flannery indicated that inner shell 
ionization and double excitation is negligible. 

The oscillator strengths calculated by Cooper's method 
were verified by means of the polarizability equation and 



24. McCann, K. J. and M. R. Flannery, "Photoionization 
of Metastable Rare-Gas Atoms," Applied Physics Letters , 
31:9:599, November 1977. 
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the calculations of Cohen and Schnieder. 



25 



They use: 



a = 



2r-2 ^ f 
e h ^ nm 



111 



^ (E -E 
m ^ n 



VI-1 



with the results listed in Table 7 . 



Table 6. Transition Array Oscillator Strengths for Ne*(3s) 



Transition 


f 


E (Ry) 
nm^ ^ ' 


Process 


3s-^2p 


-0. 0482 


1.230 


Emission 


3p 


1.080 


0.146 


Absorption 


4p 


0.0106 


0.259 


Absorption 


5p 


0.00129 


0.0302 


Absorption 



By comparison, summing over the oscillator strength of 
Table 6 and adjusting for one electron value gives: 

a = 32.492 ~ 32.5 (xlO"^'^cm^ ) . 

This agrees well with the calculations of Table 7 and demon- 
strates the accuracy attainable with the Etchechury model. 





Table 7. 


Polarizability of Ne* 


State 


(Ne*) (xl0~^^cm^) 




28.2 


Pi 


29.2 


"0 


29.2 


^Pl 


33.4 



25. Cohen, S. , and B. Schneider, "Ground and Excited 
States of Ne 2 and Ne+ 2 ," Journal of Chemical Physics , 
61:8, p. 3238, 15 October 1974. 
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The contribution to the index of refraction of the discrete 

Z f 

. The summation term indicates 

E 2-e2 

that any term having a transition energy near the laser 

cavity resonance radiation frequency (.258 Ry, 3520 A) will have 
a significant contribution to the index of refraction. From 
Table 6, it is apparent that the energy of the 3s-»-4p transition 
(E^ = .259 Ry) is near the cavity resonance energy, and re- 
quires consideration of the individual lines that comprese the 
transition array. This implies that the sum over the discrete 
transition terms in equation (11-31) is as shown in Table 8. 



Table 8. Contribution of Discrete Transitions to n-1 



Transition 




Terms in Series 


3s->-2p 


Use 


array term (negative for emission) 


3p 


Use 


array term 


4p 


Use 


summation over 30 lines 


5p 


Use 


array term 


• • • . 


Above 5p included in continuum integral 



In the transition array, there are forty possible transi- 
tions between the four 3s levels (IS 2 , Is^, Is^, in the 

Paschen notation) and the ten 4p levels (3p^, where i=l,10). 
Because of selection rules on J, ten of the lines are eliminated, 
leaving thirty lines for consideration. The relative line 
strengths for the transitions can be determined from Table 4 of 
Semenov and Strugach. The relative line strengths are 



26. Semenov, R. I., and B. A. Strugach, "Determination of 
Coupling Coefficients for the np^n's and np^n'p Configurations 
from Experimental Data," Optics and Spectroscopy , 24:, p. 260, 
April 1968. 
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indicated as the upper number and the corresponding oscillator 
strength is the lower number in Table 7. The levels of the 4p 
configuration are arrayed in order of decreasing energy across 
the top. The 3pg level is at lower energy than the 3pg. 



Table 9. Relative line strengths (upper number) and oscillator 
strengths (lower number) for 3s-^-4p transition. 



K 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 








CD 


00 




% 


3P3 


^5 


3P2 


3P4 


3p^ 


J 




.056 


0 


0 


0.139 


.833 


.111 


.667 


.333 


.833 


.333 






.0015 


j 




0.0036 


.0216 


.0029 


.0173 


.0066 


.0216 


.0066 


4 




.056 


0 


0 


.278 


0 


0 


.333 


.333 


0 


0 




2^ 


.0044 


J 


J 


.0216 


J 


J 


.0259 


.0259 


J 


J 


3 




.250 


0 


1.375 


.625 


.083 


.222 


0 


.125 


.208 


.111 


























2 




.0065 


J 


.0356 


.0162 


.0021 


.0057 




.0032 


.0054 


.0029 






.694 


2.333 


.292 


.097 


.750 


0 


0 


.208 


.625 


0 


1 


.0108 


.0361 


.0046 


.0016 


.0116 


J 




.0032 


.0097 


J 





Note: Paschen notation is used to identify levels. 

The indices K, J refer to the transitions listed in Fig. 7. 



In order to clarify the contribution of the 3s->4p transi- 
tion to the index of refraction, we consider the summation term 
in equation (11-31) in terms of cavity wavelength: 



4'(E) = C 2*°” 



nra 



Ne^ 

2 Time 



S nm o 

1^x2 

nm 



nm 



VI-2 



Figure 6 is a plot of function VI-2, using values for f 



nm 



from Table 9. Line width and damping are not considered in 
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these calculations. 



Without these considerations, the contribu- 



tions to n approaches infinity at line resonance. 

Examination of Figure 7 indicates anomalous dispersion at 
351.2, 351.6, and 351.2. The 2s^->-3pg transition is the source 
at 351.6, the 2s2^3p^ transition at 352.1, and the 2sg-^3p^Q 
transition at 352.1. 

As before, the discrete cross sections above 5p level become 
merged and indistinct, and the continuum cross section data is 
extrapolated into this region. A polynomial was fitted to the 
calculations of McCann and Flannery to generate a function for 
the integrand of equation 11-28 of the form of the previous 
section. The data and the polynomial approximation used in the 
region near the peak cross section are indicated in Fig. 4. 
Ionization threshold for Ne*(3s) is 0.365 Ry and the lower limit 
of integration is extended to 0.360 Ry, determined by the energy 
of the 5p level. The form of the integral is: 



The integration extends only to 3.3 Ry, due to the fact the 
results of McCann and Flannery terminate at that energy. The 
question arises as to the extent of error introduced by neg- 
lecting the integral from 3.3 Ry to infinity. Previous calcu- 
lations demonstrate the order of magnitude of the error intro- 
duced by termination of the integral nearer the ionization 
level. Since data are not available for the cross section 
beyond 3.3 Ry level, an extrapolation of the relationship 



00 




. 360 



th 



.395 
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0.03 
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Figure 6. Contribution of 3s-^4p Transition in Neon* to 6 ' as a Function of Wavelength 
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Figure 7. Sample Computer Output for the Calculation of the 
Contribution to the Index of Refraction of the 
3s->4p Transition. 
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Figure 8. Sample Computer Output for the Calculation of the 
Contribution to the Index of Refraction of the 
3s^4p Transition. 
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found in the calculations for unexcited neon was used. This 
implies that the cross section curves of Neon I and Ne*(3s) 
have similar shapes characterized by a peak near the ioniza- 
tion threshold and a sloping off to an asymptotic minimiim value. 
Intuitively, this approach appears correct if we compare the 
shape of the curves near the ionization threshold. A more 
rigorous discussion was not attempted. 

From Section V, the contribution to the continuum integral 
of the region near the peak is .73 of the total contribution. 
Modifying equation (11-31) by including the weighting factor 



6(E) = C 



Z 



f 

nm 





3.3 



+ BR 



/ 

.360 



g(E' )dE» 
E' ^-E^ 



VI-3 



where B 
index of 



(l/.73)= 1.37. 
refraction from 



Calculation of the variation of the 
unity, using equation (VI-3) yields 
6(E) = -7.79 X 10"^°. 
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VII. COMPARISON OF RESULTS OF CALCULATION 



The phase shifts across the cavity aperture per unit length 
of cavity may be obtained using the relation 

^ VII-1 

where L and X are expressed in consistent units. A lag in 
phase relative to vacuum is treated as a negative phase angle. 
Evaluating equation (VH-l) at 3520 A and introducing the value 
of 5(E) calculated in the previous section yields; 

^ = 2.213x10"^ m"^ 

Accuracy of the calculation is limited by the available cross 

section data beyond the ionization threshold. Improvement in 

the calculation would result from extension of the cross section 

calculations used in Section VI to an as 3 rmptotic mindmum value, 

as has been done in the case of Neon I. 

27 

Hunter , in a separate approach to the problem, calculates 

—5 —1 

the phase shift due to Ne* to be 2x10 cm 

The model developed in this paper demonstrates good agree- 
ment with independent research and is applicable to any atomic 
species present in the XeF cavity. The accuracy is dependent 
on (a) the availability of cross section data beyond the ioni- 
zation threshold of the species under consideration, and (b) 
estimates of the species population in the cavity. Current 



27. Hunter, R. 0., "Draft on Transient Index of Refraction 
in HgCl and Xe? Lasers," unpublished report, Western Research 
Corporation, Los Angeles, CA, p. 17, April 1978. 
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research into development of a populations model 
conducted by LT Lonnie Cole as a continuation of 



is being 
this research. 
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